REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  V A  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

01-04-2009  Journal  Article 


4.  TITLE  AND  SUBTITLE 

DOA  estimation  and  multi-user  interference  in  a  two-radar  system 


6.  AUTHOR(S) 

Maria  Greco*,  Fulvio  Gini*,  Alfonso  Farina**,  and  Muralidhar  Rangaswamy 


3.  DATES  COVERED  (From  -  To) 

March  2007  -  March  2008 


5a.  CONTRACT  NUMBER 

N/A 

5b.  GRANT  NUMBER 

N/A 

5c.  PROGRAM  ELEMENT  NUMBER 

612311 

5d.  PROJECT  NUMBER 

2311 

5e.  TASK  NUMBER 

HE 

5f.  WORK  UNIT  NUMBER 

02 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Electromagnetic  Scattering  Branch  *Univ.  of  I 

Air  Force  Research  Laboratory,  Sensors  Directorate  **  Selex-S 
80  Scott  Drive  Rome, 

Hanscom  AFB,  MA  01731-2909 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Electromagnetics  Technology  Division  SC:  437490 

Sensors  Directorate 

Air  Force  Research  Laboratory 

80  Scott  Drive  Hanscom  AFB,  MA  01731-2909 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Distribution  A:  Approved  for  UNLIMITED  DISTRIBUTION 


*Univ.  of  Pisa,  Italy 
**  Selex-Sistemi  Integrati 
Rome,  Italy 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR  S  ACRONYM(S) 

AFRL-RY-HS 


11.  SPONSOR/MONITOR  S  REPORT 
NUMBER(S) 

AFRL-RY -HS-TP-2009-0006 


13.  SUPPLEMENTARY  NOTES 

Research  support  by  Air  Force  Office  of  Scientific  Research.  Cleared  for  public  release  approved  by  88  ABW/PA  Public  Affairs 
Office,  WPAFB,  OH:  RY-08-0423 


14.  ABSTRACT 

In  this  paper, we  investigate  the  impact  of  the  presence  of  interfering  radar  on  the  target  direction  of  arrival  (DOA)  estimation 
performed  by  the  reference  radar. The  analyzed  estimators  are  the  pseudo-monopulse  and  the  maximum  likelihood  techniques. The 
importance  of  the  use  of  codes  in  multi-user  radar  system  is  highlighted  in  a  simple  scenario  of  two  radars  by  calculating  the 
root-mean  square  error  of  the  estimators  in  different  operational  conditions  and  comparing  them  with  the  Cramer- Rao  lower  bounds. 


15.  SUBJECT  TERMS 

Frequency-hop-coded  signals,  auto  and  cross-ambiguity  function,  monopulse  radar,  DOA  estimation,  maximum  likelihood, 
Cramer- Rao  lower  bounds,  radar  networks 


|  16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

U 

U 

U 

uu 

11 

Muralidhar  Rangaswamy 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

781-377-3446 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


Signal  Processing  89  (2009)  355-364 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Signal  Processing 

journal  homepage:  www.elsevier.com/locate/sigpro 


SIGNAL 

PROCESSING 


DOA  estimation  and  multi-user  interference  in  a  two-radar  system^ 

Maria  Greco3*,  Fulvio  Gini3,  Alfonso  Farina13'* 1,  Muralidhar  Rangaswamy c'2 

a  Dip.  di  Ingegneria  dell’Informazione,  University  of  Pisa,  via  G.  Caruso,  14-56122  Pisa,  Italy 
b  SELEX-Sistemi  Integrati,  via  Tiburtina  Km.  12.400-00131  Rome,  Italy 

c  Air  Force  Research  Laboratory  Sensors  Directorate,  80  Scott  Dr.  Hanscom  Air  Force  Base,  MA  01731-2909,  USA 


ABSTRACT 


In  this  paper,  we  investigate  the  impact  of  the  presence  of  interfering  radar  on  the  target 
direction  of  arrival  (DOA)  estimation  performed  by  the  reference  radar.  The  analyzed 
estimators  are  the  pseudo-monopulse  and  the  maximum  likelihood  techniques.  The 
importance  of  the  use  of  codes  in  multi-user  radar  system  is  highlighted  in  a  simple 
scenario  of  two  radars  by  calculating  the  root-mean  square  error  of  the  estimators  in 
different  operational  conditions  and  comparing  them  with  the  Cramer-Rao  lower  bounds. 
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1.  Introduction 

For  many  years,  conventional  radars  transmitted, 
received,  and  processed  the  same  waveform  on  every 
pulse  or  burst  within  a  coherent  processing  interval, 
independently  of  the  environment.  In  the  1970s,  adaptive 
processing  began  to  be  developed.  For  the  first  time,  the 
processing  of  received  signals  changed  depending  on  the 
environment  (noise,  clutter,  and  interferences).  Radars 
began  to  be  more  flexible  on  receive. 

Now,  modern  radar  systems  have  considerable  flex¬ 
ibility  in  their  modes  of  operation,  both  on  receive  and 
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transmit.  In  particular,  it  is  possible  to  modify  the 
waveform  on  a  pulse-to-pulse  basis,  and  electronically 
steered  phased  arrays  can  quickly  point  the  radar  beam  in 
any  feasible  direction.  Such  flexibility  calls  for  new 
methods  of  designing  and  scheduling  the  waveforms  to 
optimize  the  radar  performance.  Then,  an  agile  and 
diverse  waveform  radar  system  should  be  able  to  change 
on  the  fly  the  transmitted  waveform  based  on  the 
information  estimated  or  a  priori  known  on  the  environ¬ 
ment,  on  the  targets  and  the  jammers  [1]. 

Moreover,  in  a  radar  network  each  sensor  should  also 
be  able  to  operate  and  perform  its  task  without  negatively 
interfering  with  the  other  sensors  and,  possibly,  to 
improve  the  performance  of  the  whole  network.  Then, 
the  waveforms  used  by  the  radars  in  a  complex  network 
should  be  designed  and  changed  on  the  fly,  based  on  the 
clutter,  target  and  interference  echoes3;  they  should 
guarantee  good  target  detection  and  parameter  estimation 


3  The  choice  of  the  transmitted  waveform  depends  not  only  on  the 
desired  delay  and  Doppler  resolution,  but  also  on  where  the  clutter  or 
competing  targets  are  located  in  the  delay-Doppler  plane  [5]. 
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in  different  scenarios  and  should  allow  an  optimal  access 
to  the  same  transmit  channel. 

To  perform  good  target  detection  and  estimation, 
frequency  hop  pulse  train  signals  are  often  used  in  high 
resolution  radar  systems.  These  signals  are  characterized 
by  an  auto-ambiguity  function  (AAF)  that  exhibits  a 
narrow  thumb  tack  shape  with  low  sidelobes.  In  contrast, 
in  application  like  multi-access  communications,  atten¬ 
tion  is  paid  in  designing  a  sequence  of  frequency-hopped 
coded  waveforms  with  small  cross-correlation  functions. 
In  multi-user  radar  system  scenarios  both  objectives  are 
desirable.  Unfortunately,  there  is  a  tradeoff  between  these 
objectives.  Frequency  hop  pulse  trains  based  on  Costas 
codes,  for  instance,  are  known  to  have  almost  ideal  AAF 
but  not  very  good  cross-ambiguity  properties  [2].  On  the 
contrary,  frequency-coded  signals  based  on  linear  con¬ 
gruences  [3]  have  ideal  cross-  but  unattractive  auto¬ 
ambiguity  properties.  Some  attempt  to  design  multiple 
access  frequency  hop  codes  with  good  cross-ambiguity 
function  (CAF)  has  been  done,  for  instance,  in  [4]  where 
the  frequency  hop  patterns  were  constructed  upon  an 
extension  of  the  theory  of  quadratic  congruences. 

The  scenario  analyzed  in  this  paper  is  composed  of  two 
radars  transmitting  in  the  same  band,  which  can  illumi¬ 
nate  the  same  area  looking  for  the  same  target;  they  can 
use  either  the  same  or  different  codes.  We  rewrite  the 
observed  signal  as  a  function  of  cross-  and  auto-ambiguity 
function  values  and  target  parameters,  and  we  investigate 
the  impact  of  the  presence  of  the  transmitted  signal  of  the 
second  radar  (the  interfering  radar)  on  the  first  one  (the 
reference  radar)  in  the  estimation  of  the  target  direction  of 
arrival  (DOA).  We  analyzed  two  DOA  estimators,  the 
pseudo-monopulse  (PM)  and  the  maximum  likelihood 
(ML)  estimator  and  we  compared  their  performances  with 


the  Cramer- Rao  lower  bounds  (CRLB).  We  focused  our 
attention  only  on  the  influence  of  the  interfering  radar  on 
the  DOA  estimation,  so  in  our  scenario  we  did  not 
simulate  jammers  or  correlated  clutter. 

The  paper  is  organized  as  follows.  In  Section  2,  the 
received  signal  model  is  explained  and  the  PM  and  the  ML 
techniques  for  DOA  estimation  are  summarized.  In  Section 
3,  the  results  of  the  analysis  are  described.  In  Section  4, 
the  CRLB  for  the  problem  at  hand  is  derived  and  compared 
with  the  root-mean  square  error  (RMSE)  of  the  PM  and  ML 
estimators.  In  Section  5,  some  conclusions  on  the  use  of 
codes  are  drawn. 


2.  Problem  statement 


The  scenario  is  pictorially  drawn  in  Fig.  1.  Two  radars 
scan  the  same  area  transmitting  frequency-coded  bursts 
of  pulses  and  listening  to  the  echo  of  a  possible  target.  The 
complex  envelope  of  the  transmitted  unitary  power  signal 
is  given  by 

-1  M 

u®  =  ^o=^Um(t-(m-1)Tc)’  (!) 

v  ™Tc  m=1 

where 


f  exp0'27i/mt)  O^t^Tc, 

1  0  elsewhere, 


(2) 


M  is  the  number  of  subpulses  for  each  transmitted  pulse 
of  time  duration  Tiy  tc  =  TilM  is  the  duration  of  each 
subpulse  and  i  is  the  sequence  of  frequencies 

related  to  the  code  used  by  the  radar.  For  Costas  arrays,  for 
instance,  fm  =  dm\ tc,  where  dm  belongs  to  the  sequence 
dM  =  .  .  -  which  is  a  permutation  of  the 


Fig.  1.  Radar  scenario. 
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integer  numbers  Jm={ 0,1,...,M-1}.  The  choice  of 
{/mWl  i,  characterizing  the  frequency  code  is  critical 
and  of  paramount  importance  in  defining  the  properties  of 
the  auto-  and  cross-ambiguity  functions  of  the  trans¬ 
mitted  signal.  A  detailed  description  of  this  topic  can  be 
found  in  [5]. 

As  known,  the  AAF  represents  the  time  response  of  a 
filter  matched  to  a  given  finite  energy  signal  when  the 
signal  is  received  with  a  delay  z  and  a  Doppler  shift  v 
relative  to  the  nominal  zeros  values  expected  by  the  filter 
[5].  Then,  the  AAF  definition  is 


|A(t,v)| 


u(t)u*(t  +  t)  exp(j27ivt)dt 


(3) 


where  u(t)  is  the  complex  envelope  of  the  signal.  The  CAF 
between  two  signals  iq(t)  and  u2(t)  is  similarly  defined  as4 


|C(T,V)| 


u.\  (t)u5(t  +  t)  exp(j27ivt)dt 


(4) 


2.1.  Signal  model 


Fig.  2.  Receiver  scheme. 

transmitted  signal  u(n//c),  we  can  state  that  A  is  the  value 
of  the  complex  AAF  of  the  signal  for  some  delay  z  and 
Doppler  target  shift /D;  then  A  =  A(t/d).  This  is  the  case  of 
the  signal  backscattered  by  a  target  illuminated  by  the 
reference  radar.  If,  on  the  contrary,  the  signal  received  by 
the  reference  radar  is  the  signal  transmitted  by  the 
interfering  radar,  A  is  the  value  of  the  complex  CAF;  that 
is  A  =  C(tr,/Dr).  If  both  signals  are  present,  the  received 
signal  is  given  by 


In  the  reference  radar,  the  received  signal  is  first  down¬ 
shifted  to  an  intermediate  frequency  IF  and  amplified.  The 
IF  signal  is  then  processed  as  in  the  scheme  of  Fig.  2, 
where  LO  is  a  local  oscillator  tuned  on  the  IF  frequency, 
LPF  is  a  low  pass  filter  and  A/D  is  an  analog-to-digital 
device. 

Before  the  digitalization,  the  inphase  (I)  and  (Q) 
quadrature  components  of  the  target  signal  are: 

Xi(t)  =  a  cos  ^271  +  2nfDt  -  (p^j,  (5) 

Xq(£)  =  a  sin  ^2ti t  +  27i/Dt  —  ,  (6) 

where  aexp(jcp)  is  the  complex  amplitude  of  the  target 
that  depends  on  the  radar-cross  section  (RCS)  and  on 
the  antenna  gain,  fD  is  the  Doppler  frequency,  and 
d(t)  =  J2m^Todm-  rect((t-mTc/2)Tc)  is  the  frequency  code. 

After  the  digitalization  at  a  sampling  frequency  /c,  the 
complex  envelope  of  the  received  target  signal  is  given  by 

x(n//c)  =  ae^  exp  ^ 2n  n  +  271  J^n^)  •  (7) 

In  the  correlator  of  Fig.  2,  the  sequence  of  x(n//c)  is 
correlated  with  the  sequence  of  samples  of  the  trans¬ 
mitted  signal.  Then  the  output  signal  is  given  by5 

Ns- 1 

x  =  ae^  ]X  u(n//c)x*(n//c)  =  ae^A.  (8) 

n= 0 

It  is  easy  to  compare  the  term  A  in  Eq.  (8)  with  Eqs.  (3) 
and  (4).  If  Ns  is  large  we  can  consider  the  sum  as  a  good 
approximation  of  the  integral.  Moreover,  if  the  received 
signal  x(n//c)  is  a  delayed  and  frequency  shifted  copy  of  the 


4  The  calculation  of  Eqs.  (4)  and  (5)  is  not  always  easy.  For  some 
frequency  code  it  is  useful  to  resort  to  the  placement  and  difference 
matrices,  as  explained  in  [3,14]. 

5  We  did  not  use  any  window  in  the  calculation  of  Eq.  (8). 


y  =  aei<M(T,/D)  +  be>0C(TR,/DE)  +  d,  (9) 

where  be*9  is  the  complex  amplitude  of  the  signal  relating 
to  the  second  radar  and  d  is  the  unavoidable  contribution 
of  the  disturbance  (thermal  noise  plus  clutter). 

2.2.  Pseudo-monopulse  estimator 

In  the  analyzed  scenario  the  reference  radar  is 
supposed  to  estimate  the  DOA  of  the  target  by  using  a 
pseudo-monopulse  technique  or  the  maximum  likelihood 
estimator  [6]. 

In  a  typical  phased  array  radar,  a  single  beam  is  formed 
on  transmission  and  two  or  more  beams  are  formed  on 
reception.  We  assume  here  that  the  system  is  a  linear 
array  radar  that  estimates  the  target  DOA  by  using  the 
sum  channel  XI  on  transmission  and  two  matched 
channels,  the  sum  X  and  the  difference  A  on  reception. 
The  two  channels,  or  antenna  patterns,  are  defined  as  the 
complex  amplitude  profiles  versus  target  azimuth  angle 
0T.  The  sum  channel  pattern  is  denoted  by  /^(0)  and  the 
difference  pattern  by  these  patterns  are  chosen  here 
as  in  [6],  the  -3dB  beamwidth  is  0B  =  3°.  The  PM 
technique  used  here  is  very  similar  to  the  classical 
monopulse  method,  based  on  the  ratio  of  the  A  and  X 
channel  outputs  yA  and  y^.  It  uses  N  pulses  transmitted  by 
the  radar  in  the  time-on-target  (ToT),  while  the  antenna 
beams  are  electronically  or  mechanically  steered  with 
constant  angular  velocity  coR  rad/s.  The  number  N  of 
pulses  between  the  one-way  -3dB  points  is  given  by 
N  =  ObI(corT),  where  T  =  1/PRF  is  the  radar  pulse  repeti¬ 
tion  time  and  PRF  is  the  corresponding  pulse  repetition 
frequency.  The  antenna  introduces  an  amplitude  modula¬ 
tion  on  the  target  signal,  in  both  the  channels,  that 
depends  on  the  target  azimuth  position  and  on  the 
instantaneous  boresight  of  the  array. 

Based  on  Eq.  (9),  it  is  possible  to  write  the  expression  of 
the  signal  received  on  the  sum  and  on  the  difference 
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channels.  Let  us  suppose  that  the  Doppler  frequencies  of 
the  target  and  of  the  interference  radar  are  the  same  with 
respect  to  the  reference  radar.  In  this  case,  the  reference 
radar  cannot  distinguish  between  the  two  signals  by 
means  of  Doppler  processing.  Without  lack  of  generality, 
we  assume  that  both  Doppler  frequencies  are  null;  then 
the  received  signals  are  given  by 

y£(n)  =  aA/|(0T,n)  +  /?C/2(0,,n)  +  d£(n),  (10a) 


is  given  by 


Oml  =  argmax 

9 


lyHM  ag(0)i2 
gH(6>)M-1g(6»)' 


(11) 


In  the  next  paragraph  we  investigate  the  impact  of  the 
presence  of  the  interfering  radar  on  the  DOA  RMSE  when 
the  estimation  is  performed  by  the  PM  and  the  ML 
techniques. 


yA(n)  =  y.AfL(0T,  rt)fA(9T,  n)  +  fiCfA(Ou  n)  +  dA(n),  (10b) 

where  n  =  0,1,..., N-l.  With  respect  to  Eq.  (9),  the 
dependence  of  the  complex  amplitudes  of  target  and 
interfering  signal  on  the  antenna  patterns  has  been 
explicitly  indicated.  The  first  term  in  and  yA  is 
due  to  the  target  signal;  it  depends  on  the  target 
DOA  0T  through  /|(0T)  in  sum  channel  and  through 
fz{0T)fA(0T)  in  the  difference  channel.  This  is  due  to  the 
two-way  antenna  gain.  The  second  term  is  due  to  the 
signal  transmitted  by  the  second  radar  which  has  a  DOA  6{ 
and  it  depends  on  the  one-way  gain  of  the  antenna 
pattern  Cfe(di)  and/^di)).  dz  and  dA  are  the  noises  on  the 
two  channels. 

The  signal  processor  forms  the  monopulse  ratio 
defined  by  r(n)  =  Re{yA(n)ly^(n)}  for  each  pulse  where 
Re{}  denotes  the  real  part.  In  absence  of  disturbance  and 
in  the  presence  of  only  one  target,  the  monopulse  ratio 
reduces  to  r(n)  =  Re{fA(0T,n)lfe{0T,n)}  from  which  the 
angular  location  of  the  target  for  each  pulse  can  be 
determined.  Finally,  we  obtain  0TG  =  Yln=oOjG/^,  where 
Ojq  is  the  target  DOA  estimate  for  each  pulse. 

We  treat  the  situation  in  which  the  source  can  be 
considered  static  with  respect  to  the  radar  during  the  ToT, 
that  is,  during  the  recording  of  the  N  pulses.  In  this  case, 
due  to  the  scanning  movement  of  the  antenna  only,  the 
proposed  estimator  is  biased  and  the  bias  can  be 
calculated.  It  is  easy  to  verify  that  b  =  E{6TG  -  0jG}  = 
(N  —  1)0r/(2N)  then  a  DOA  unbiased  estimator  is  0TG  = 
J2nX /N-b  [6]. 


2.3.  Maximum  likelihood  estimator 

If  we  rewrite  Eq.  (10)  in  vector  notation,  the  data  model 
is  given  by  y  =  aig(0T)+aif(0i)+d,  where  aT  =  aA,  and 
(Xi  =  fiC,  g (0T)  =  [gsgjT]T  is  the  2N  x  1  steering  vector  of 
the  target,  [g^]„  =fi(0T,n)  and  [gA]n  =  M0T,n)fA{0T,n). 
f(0T)  =  [fzf/]T  is  the  steering  vector  of  the  interference, 
with  [fz]n  =/2(0T,n)  and  [fA]n  =/zi(0T,n).  The  2Nxl  dis¬ 
turbance  vector  d  is  given  by  d  =  [d^dAT]T.  The  ML 
estimator  of  0T,6  when  the  observed  vector  does  not 
contain  the  interference  and  the  vector  d  is  modeled 
as  a  complex  Gaussian  vector  with  a  covariance  matrix  M, 


6  As  shown  in  [6],  Eq.  (11)  represents  the  ML  estimator  of  0T  when 
the  amplitude  aT  is  modeled  as  a  deterministic  unknown  value.  When  aT 
is  a  random  variable,  the  ML  estimator  has  not  a  closed  form.  We  adopt 
here  the  ML  estimator  for  deterministic  aT  even  if  in  our  data  model  the 
complex  amplitude  is  random;  with  an  abuse  of  terminology  we 
continue  to  refer  to  (11)  as  the  ML  estimator. 


3.  Simulation  results 

To  evaluate  the  impact  of  the  presence  of  the  interfer¬ 
ing  radar,  the  RMSE,  the  variance  and  the  bias  of  the  DOA 
estimator  has  been  derived  by  running  104  Monte  Carlo 
simulations.  The  disturbance  is  modeled  as  a  complex 
white  Gaussian  process,  in  short  notation  d~CN(0,Odl), 
where  I  the  2N-dimensional  identity  matrix,  the  target 
and  interference  signal  amplitudes  aT  and  oq  are  modeled 
first  as:  (i)  complex  Gaussian  independent  random 
variables,  in  short  aT~CN(0,cr2|A|2)  and  ar~CN(0,cr2|C|2); 
and  then  as:  (ii)  deterministic  unknown  parameters.  In 
the  first  case,  the  signal-to-noise  (SNR)  power  ratio  is 
defined  as  SNR  =  |A|2  •  o\\o\  and  the  signal-to-interfer- 
ence  ratio  as  SIR  =  cr2/cr2;  in  the  second  case,  SNR  = 
|aT|2  •  /< 7d  and  SIR  =  |aT/ai|2. 

The  results  of  our  analysis  are  shown  in  Figs.  3-7  for 
the  PM  estimator  and  in  Figs.  8-14  for  the  ML  estimator. 
We  set  SNR  =  20dB,7  SIR  =  OdB  and  N  =  16;  the  target 
DOA  has  been  set  equal  to  1.5°  (0T  =  1.5°).  We  investigated 
also  other  target  DOA  positions  in  the  main  beam  of  the 
antenna.  The  behavior  of  the  estimators  is  similar,  and 
then  we  show  here  the  results  for  one  position  only. 

Due  to  the  high  value  of  SNR,  the  performance  of  the 
PM  and  ML  estimators  are  mainly  affected  by  the  presence 
of  the  interfering  radar;  they  are  the  functions  of  the  ratio 
|  C/A | .  The  value  of  C  and  A  depends  on  the  code  used  by 
both  radars  and  on  the  synchronization  between  transmit 
and  receive  of  the  reference  radar  and  between  reference 
and  interfering  radar. 

If  the  reference  radar  is  synchronized  in  reception  and 
transmission  and  the  receiver  is  tuned  on  the  Doppler  of 
the  target,  t  =  0  and  fD  =  0,  then  A  is  the  energy  of  the 
transmitted  pulse.  In  our  analysis  we  considered  A  =  1 
(unit  energy  pulse)  and  0^|C/A|^1.  The  worst  case  is 
when  |C/A|  =  1;  this  value  characterizes  two  synchronized 
radars  using  the  same  code.  The  best  case  is  when  |C/ 
A|  =  0;  this  value  characterizes  the  case  of  synchronized 
radars  using  two  orthogonal  codes. 

In  Fig.  3,  we  report  the  RMSE(0T)  of  the  PM  estimator  as 
a  function  of  the  interference  DOA,  with  aT~CN(0,<72|A|2) 
and  ai~CN(0,cr2|C|2).  In  Fig.  4,  we  report  the  bias  of  the  PM 
estimator  in  the  same  scenario.  The  bias  of  the  estimator 
increases  with  |C/A|  and  the  RMSE  gets  values  close  to  the 
40%  of  the  beamwidth.  The  behavior  of  both  RMSE  and 
bias  are  almost  symmetric  with  respect  to  the  position  of 
the  target  0T.  When  =  0T,  the  bias  due  to  the  presence  of 


7  With  this  value  of  SNR  the  probability  of  detection  of  the  filter 
matched  to  the  transmitted  signal  is  almost  unitary,  even  for  PFA  =  10-6. 
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Fig.  3.  RMSE  vs.  interference  DOA,  0T=1.5°,  SNR  =  20dB,  SIR  =  0dB, 
N  =  16,  random  a  and  /?,  PM  estimator. 


0j  (degrees) 

Fig.  4.  Bias  vs.  interference  DOA,  0T=  1.5°,  SNR  =  20dB,  SIR  =  0dB, 
N  =  16,  random  a  and  /?,  PM  estimator. 
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Fig.  6.  Bias  vs.  interference  DOA,  0T  —  1.5°,  SNR  =  20dB,  SIR  =  0dB, 
N  =  16,  deterministic  gc  and  /?,  PM  estimator. 
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Fig.  7.  RMSE  vs.  |C/A|,  0T  =  1.5°,  SNR  =  20dB,  N=16,  0I£D(-7O,7O), 
random  a  and  /?,  PM  estimator. 
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Fig.  5.  RMSE  vs.  interference  DOA,  0T=l-5°,  SNR  =  20dB,  SIR  =  0dB, 
N  =  16,  deterministic  a  and  /?,  PM  estimator. 
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Fig.  8.  RMSE  vs.  |C/A|,  0T  =  1.5°,  SNR  =  20dB,  N=16,  0iel/(- 7°, 7°), 
deterministic  a  and  /?,  PM  estimator. 
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Fig.  9.  RMSE  vs.  interference  DOA,  0T=  1.5°,  SNR  =  20dB,  SIR  =  OdB,  Fig.  12.  Bias  vs.  interference  DOA,  0T=  1.5°,  SNR  =  20dB,  SIR  =  OdB, 

N  =  16,  random  a  and  p,  ML  estimator.  N  =  16,  deterministic  a  and  p,  ML  estimator. 


0j (degrees) 


Fig.  10.  Bias  vs.  interference  DOA,  0T=  1.5°,  SNR  =  20dB,  SIR  =  0dB,  Fig.  13.  RMSE  vs.  |C/A|,  0T  =  \.5°,  SNR  =  20dB,  N=16,  0iel/(-7°,7°), 

N  =  16,  random  a  and  p,  ML  estimator.  random  a  and  p,  ML  estimator. 
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Fig.  11.  RMSE  vs.  interference  DOA,  0T  =  1.5°,  SNR  =  20  dB,  SIR  =  OdB, 
N  —  16,  deterministic  a  and  p,  ML  estimator. 


Fig.  14.  RMSE  vs.  |C/A|,  0T  =  1.5°,  SNR  =  20dB,  N=16,  0ieD(-7°,7°), 
deterministic  a  and  p,  ML  estimator. 
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the  interfering  radar  is  null  and  the  RMSE(0T)  for  |C/A|  =  1 
is  lower  than  that  for  |C/A|  =  0. 

Similar  remarks  can  be  drawn  from  Figs.  5  and  6, 
where  we  report  RMSE(0T)  and  bias  for  deterministic 
unknown  aT  and  ai.  For  small  |C/A|,  the  RMSE(0T)  of  the  PM 
estimator  is  lower  when  the  target  and  interference 
amplitudes  are  deterministic  than  when  they  are  complex 
Gaussian  random  variables.  Conversely,  the  RMSE(0T) 
is  more  sensitive  to  high  values  of  |C/A|  =  0  when 
the  amplitudes  are  deterministic  than  when  they  are 
random. 

In  Figs.  7  and  8,  RMSE(0T)  is  shown  as  a  function  of  the 
ratio  |C/A|  for  different  value  of  the  SIR,  in  the  case  of 
random  and  deterministic  amplitudes,  respectively.  The 
interference  DOA  0X  has  been  generated  as  a  random 
variable  uniformly  distributed  on  the  range  [-7°,-7°].8 
The  impact  of  the  interference  is  apparent  in  both 
figures. 

Figs.  9-14  report  the  companion  results  for  the  ML 
estimator.  The  behavior  is  very  similar.  In  the  case  of 
deterministic  amplitudes  the  ML  estimator  outperforms 
the  PM  estimator;  the  contrary  in  the  random  amplitude 
scenario.  It  is  worth  reminding  that  in  the  case  of  random 
amplitude,  what  in  this  paper  is  termed  “ML  estimator”  is 
not  the  true  ML  estimator  (for  details  see  [7]). 


4.  Cramer-Rao  lower  bounds  and 
performance  comparison 


To  complete  our  analysis  we  derived  also  the  CRLB  for 
the  two  cases  at  hand.  In  this  section  we  report  results  and 
comments.  Detailed  derivation  of  the  bounds  is  presented 
in  Appendices  A  and  B. 

When  the  observed  vector  y  is  complex  Gaussian 
distributed  with  mean  value  py  and  covariance  matrix  Cy, 
in  short  y~CN(py,Cy),  with  unknown  real  parameters 
X  =  [XbX2, •••Xn]T»  the  elements  of  the  Fisher  matrix  are 
given  by  [8] 


D(X)]«  =  Trjc;1® 


f  011?  (7)  0U, 

-2s 

\  3 Xt  y  dXj 


(X) 


(12) 


where  Tr{  • }  is  the  trace  of  a  matrix,  91e{  • }  stands  for  the 
real  part  and  ij  =  1,2 ,...,n. 


4.1.  Deterministic  amplitudes 

In  this  case  Cy  =  ail  and  py  =  otA  -g(0T)+^C- h(0j). 
Decomposing  the  complex  amplitudes  a  and  in  their 
modulus  an  phase,  that  is  a  =  ae1^  and  /3  =  he1^,  the  vector 
of  the  real  parameters  to  estimate  is  /  =  [0T  a  (p  0\  b  il/]T. 
The  36  elements  of  the  Fisher  matrix  are  reported  in 
Appendix  A. 

As  verified  in  Eq.  (A.28),  the  CRLB  of  the  target 
parameters  0T,  a,  and  cp  do  not  depend  on  the  ratio  |C/A|. 


This  means  that,  theoretically,  if  the  radar  knows  that  an 
interference  is  present,  it  is  possible  to  find  an  efficient 
unbiased  estimator  that  is  not  sensitive  to  the  code  used 
by  both  reference  and  interfering  radars.  Unfortunately, 
this  never  happens.  In  the  existing  systems,  the  radar 
simply  applies  some  techniques  for  the  estimation  of 
target  DOA  even  when  there  is  an  interference,  then  the 
presence  of  an  interfering  radar  biases  the  estimates 
increasing  the  root  mean  square  error  more  than  one 
order  of  magnitude  with  respect  to  the  CRLB,  as  shown 
in  Fig.  15.  In  Fig.  15,  we  report  the  square  root  of 
the  Cramer-Rao  lower  bound  (RCRLB)  of  0T  in  degrees, 
RCRLB(0t)  =  (1 80°  /71)  CRLB(0t),  compared  to  the 
RMSE(0T)  of  ML  and  PM  estimators  for  |C/A|  =  0  and  |C/ 
A  |  =  1,  as  a  function  of  the  interference  DOA.  The  solid  line 
with  white  squares  (labelled  RCRLB)  represents  the 
RCRLB(0t)  calculated  when  all  the  target  and  interference 
parameters  are  deterministic  and  unknown.  The  dotted 


0j  (degrees) 

Fig.  15.  RMSE  of  ML  and  PM  estimators  and  RCRLB  vs.  interference  DOA, 
0T  —  1.5°,  deterministic  a  and  /L 
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8  These  are  the  angular  positions  of  the  nulls  of  first  right  and  left 
sidelobes  of  the  sum  channel. 


Fig.  16.  RMSE  of  ML  and  PM  estimators  and  RCRLB  vs.  interference  DOA, 
6t  =  1.5°,  random  a  and  /L 
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line  with  white  circles  (labelled  RCRLBwi)  reports 
the  RCRLB(dT)  calculated  when  the  interference  is  not 
present  (Eq.  (A.29)).  It  should  not  seem  to  be  an 
absurd  that  for  some  values  of  the  RMSE(dT)  is  smaller 
than  the  RCRLB(dT),  because  both  ML  and  PM  estimators 
are  biased. 


4.2.  Random  amplitudes 

The  comparisons  in  the  case  of  complex  Gaussian 
distributed  amplitudes  are  summarized  in  Fig.  16.  Again, 
the  solid  line  with  white  squares  (labelled  RCRLB@|C/ 
A\  =  1)  represents  the  RCRLB(dT)  calculated  when  both 
target  and  interference  DOAs  are  unknown  and  |C/A|  =  1. 
The  dotted  line  with  white  circles  (labelled  RCRLBwi) 
reports  the  RCRLB(dT)  calculated  when  the  interference  is 
not  present  (Eq.  (B.3)).  The  two  curves  are  close,  meaning 
that  also  in  this  case,  knowing  the  presence  of  the 
interference,  it  could  be  possible  to  find  an  estimator  that 
is  not  too  sensitive  to  the  values  of  |C/A|.  Unfortunately, 
also  in  the  case  of  random  amplitude,  the  performance  of 
realistic  estimators  as  the  PM  and  ML  are  far  from  the 
CRLB  and  the  impact  of  the  codes  is  heavy. 

5.  Conclusions 

The  aim  of  this  paper  is  to  highlight  that  the  proper  use 
of  codes  is  of  paramount  importance  even  in  a  very  simple 
radar  network  formed  by  two  radars.  The  impact  of  the 
presence  of  a  radar  on  the  other  has  been  measured  in  the 
estimation  of  the  target  DOA  performed  by  the  reference 
radar.  We  considered  only  two  multi-pulses  algorithms, 
the  ML  and  the  PM,  among  all  the  DOA  estimation 
techniques  (see  for  instance  [9-12],  and  references  there¬ 
in)  not  with  the  purpose  of  finding  the  best  algorithm  in 
the  presence  of  an  interference,  but  simply  willing  to 
show  that  two  very  different  techniques  suffer  the  same 
performance  impoverishment.  We  verified  that  the  inter¬ 
ference  mainly  introduces  a  bias  in  the  DOA  estimate  that 
depends  on  the  position  of  the  interference  with  respect 
to  the  boresight  of  the  antenna  and  on  the  ratio  |C/A|,  that 
is,  on  the  choice  of  the  codes.  For  limiting  this  negative 
effect,  the  code  should  be  chosen  such  that  |C/A|  <  1,  that  is 
the  codes  used  by  the  two  radars  should  be  almost 
orthogonal.  This  is  not  always  an  easy  task.  Generally  the 
code  is  designed  for  optimizing  the  frequency  and  range 
resolutions  of  the  radar,  that  is,  for  optimally  shaping  the 
AAF.  But  often  the  best-shaped  auto-ambiguity  corre¬ 
sponds  to  a  very  poor  CAF,  that  is,  in  our  problem,  high 
value  of  |  C/A | .  Some  asymptotic  bounds  on  CAF  of  different 
classes  of  frequency  hop-coded  signals  have  been  reported 
in  [3,13-15].  For  Welch-Costas  code  of  length  N  the  bound 
is  N/N  then  |C/A|  can  easily  get  the  unit  value.  For  codes 
based  upon  linear  congruences  the  bound  is  2/N,  then  the 
maximum  value  of  |C/A|  is  0.5.  Good  results  can  be 
obtained  with  a  set  of  the  codes  based  upon  extended 
quadratic  congruences  for  which  the  asymptotic  bound 
seems  to  be  12/N.  The  adaptive  use  of  different  trans¬ 
mitted  waveforms  to  alleviate  the  impact  of  an  inter¬ 
ference  is  now  under  analysis  [16]. 


Appendix  A.  CRLB  derivation  for 
deterministic  amplitudes 


First,  let  us  define  some  useful  vectors  and  matrices 


m  = 


dg  (0) 
W  9 


m  = 


6h(d) 

W 


and 


(A.1) 


G(0)  =  g(0)gH(0)  =  g(0)gT(0),  H(0)  =  h(0)h>)  =  h(0)hT(0), 

(A.2) 


G(0)  = 


0G  (0) 
“60“ 


H(0)  = 


9H(0) 

~d0~ 


(A.3) 


where  the  n.mth  element  of  the  matrix  is  [G(0)]„  m  = 
[mUS(d)]m  +  [g(0)]„[g(0)]m.  and  [H(0)]nm  =  [h(0)]n[h(0)]m 

+[h(6>)]n[h  mm. 

In  the  case  of  deterministic  amplitudes  Eq.  (12) 
reduces  to 


U(x)]y  =  4^ 

°d 


3Xi  dXj  ]’ 


ij=  1,2,...  6. 


(A.4) 


Knowing  that  fiy  =  ae^A  •g(dT)+fiei^C-h(dI),  we  can 
calculate 


OflyOO 


'  ae^’Agdh) 

e^Ag(6>T) 

jae^AgCd-r) 

fae^Ch(di) 

e^Ch(0,) 

jt>e^Ch(0i) 


(A.5) 


Setting  now  y  =  ]C/A[,  A  =  |A|ei9'1  and  C  =  |C|ei9c  = 
y|A|ei9c,  we  obtain  for  the  Fisher  matrix  elements 


Ju  =4aWHS(0T)ll2 


(A.6) 


Jn  =  J21  =  ^a|A|2gT(0T)g(0T). 


J 13  —  J31  —  0, 


(A.7) 

(A.8) 


J 14  =  J41  =  ^2  aby\Af  cos(9c  -9A  +  il/-(p)g  (0T)h(0i), 


(A.9) 


J15  =Js  1  =  ^2ay\A\2  cos  (9C  -9a  +  iA  -  <P) gT(0T)h(0i), 

(A.10) 


Jie  =Je  1  =  -  —jObylAj2  sin(3c  -  9A  +  i// -  cp)gT(6T)h(0i), 

(A.ll) 

J22=^|A|2||g(0T)||2,  (A.1 2) 


J 23  —  J32  ~  0. 


(A.13) 
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J24  =  J42  =  ^2by\A\2  cos (9C  -  9„  + 1 // -  (p)gT(0T)h(0i), 

(A.14) 

J25  =Js 2  =  ^rlA|2  cos (9C  -  9a  +  i//  -  cp)gJ (6j)h(6[), 

(A.15) 

J26  =Je 2  =  -^by|A|2  sin(9c  -9A+\ // -  cp)gT (0T)h(0,), 

(A.16) 

J33=^aW||g(0T)||2,  (A.17) 

J34  =  J43  =  ^aby\A\2  sin (9c  -9A  +  xj/  -  (c)gT(9T)h(9I), 

ffd 

(A.18) 

J35  =  J53  =  ^2a7\A\2  sin(9c  -9/1  +  1 A  -  <p)gT(0T)h(0,), 

(A.19) 


J36  =ie3  =  ^aby|A|2  cos (9c  -9/ 1  +  lA  -  (o)gT(9T)h(0i), 


(A.20) 

J44  =  ^bV|A|2||h(0,)||2, 

(A.21) 

J45  =Js4  =  by2  |A|2hT(0,)h(0,), 

(A.22) 

J  46  =J  64  =  0* 

(A.23) 

J55=^yW||h(0,)||2> 

(A.24) 

J56  =  J  65  —  0* 

(A.25) 

;66=^bV|A|2||h(9,)||2.  (A.26) 

Based  on  Eqs.  (A.6)-(A.26),  we  can  observe  that  the 
Fisher  matrix  can  be  written  as 


J  = 


A 

yBT 


yB 

y2C 


(A.27) 


where  A,  B  and  C  are  three-dimensional  sub-matrices. 
Then,  the  inverse  of  J  can  be  expressed  as 


D  E 

(A  —  BC-1Bt)-1 

-f<A  BC  'Brl  'BC  ' 

F  G 

-i(C-BTA ^B)  :BA  1 

J+C-B’A-'Br1 

(A.28) 

It  is  worth  noting  that  the  sub-matrix  D  does  not 
depend  on  7,  then  the  CRLB  of  0T,  a ,  and  cp  do  not  depend 
on  7.  Conversely,  the  bounds  of  the  interference  para¬ 
meters  are  inversely  proportional  to  y2. 

If  the  interference  is  not  present,  the  CRLB  of  0T,  a,  and 
cp  can  be  calculated  simply  inverting  the  matrix  A.  After 


some  algebra  we  get 
CRLB(9t)  =  IT1],., 

=  _J _ llg(8T)ll2 _ 

2SNR I  lg(9r)ll2 1  lg(0T)l  I2  -  gT(0T)g(0T)gT(0T)g(0T)  ’ 

(A.29) 


Appendix  B.  CRLB  derivation  for  complex  Gaussian 
distributed  amplitudes 

In  this  case,  supposing  aeCN(0,cra2)  and  ^eC^O.cr^2), 
the  observed  vector  is  still  Gaussian  with  a  mean 
value  py  =  £{y}  =  0  and  a  covariance  matrix  Cy  = 
o2  |A  |2G(  O^+Gp2 1  C|  2H(  0/)+crdl.  The  unknown  parameter 
vector  is  %  =  [0T  6i]J  and  (12)  reduces  to 


B(X)]y  =  Tr- 


Cy'W 


8Cy(X) 

3^,' 


C/(X) 


8Cy(X)| 

0  xj  y 


i,j=  1,2. 

(B.l) 


It  is  easy  to  verify  that 


8C y 

0Xi 

and 

0Cy 

0x2 


0Cy 

00T 


fT2|A|2G(9T) 


0Cy 

00, 


=  (t«ici2h(0,). 


(B.2) 


Replacing  (B.2)  in  (B.l)  and  inverting  the  Fisher  matrix 
we  can  obtain  the  CRLB. 

If  the  interference  is  not  present,  the  Cramer-Rao 
bound  of  0T,  can  be  easily  calculated  inverting  [J]1/L.  After 
some  algebra 

CRLB(0t) 

1  1 

“  SNRir{[(I  -  (SNR  •  G(<9T)/1  +  SNR  •  gT(0T)g(0T)))  •  G(3T)]2} ' 

(B.3) 
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